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I. TURBULENCE FOR THE PHYSICIST AND 
FOR THE ENGINEER 

Sir Horace Lamb once said "I am an old man now, 
and when I die and go to Heaven there are two mat- 
ters on which I hope enlightenment. One is quantum 
electro-dynamics and the other is turbulence of fluids. 
About the former, 1 am really rather optimistic" [Q. 
Possibly Lamb's pessimism about turbulence had been 
short-sighted. There exist signs that the two issues that 
concerned Lamb are not disconnected. The connections 
were brewing for some while, and began to take clearer 
form recently. This promises renewed vigor in the in- 
tellectual pursuit for understanding this long-standing 
problem. This is not due to some outstanding devel- 
opments of new tools in the theoretical or experimental 
study of turbulence per se, but rather due to develop- 
ments in neighboring fields. The great successes of the 
theory of critical and chaotic phenomena and the popu- 
larity of nonlinear physics of classical systems attracted 
efforts that combined the strength of fields like quantum 
field theory and condensed matter physics leading to re- 
newed optimism about the solubility of the problem of 
turbulence. It seems that this area of research will have 
a renaissance of rapid growth that promises excitement 
well into the next century. 

It is possibly a happy coincidence that hydrodynamic 
turbulence is considered a problem of immense interest 
by both physicists and engineers. The physicist tends to 
appreciate phenomena that display universal characteris- 
tics; the engineer may find such characteristics irrelevant 
since they cannot be manipulated. The engineer seeks 
control, and control means a ready response to pertur- 
bations. Universal phenomena are immune to perturba- 
tions. The point is of course that "turbulence" means 
different things to different researchers. All agree that 
hydrodynamic turbulence arises in fluids that are highly 
stressed, or stirred, such that there exist significant fluid 
velocities (or winds) on the largest scales of motion. The 
engineer is typically interested in the flow characteristics 
near the boundaries of the fluid (boundary layers, air- 
plane wings, pipes, turbines etc). By understanding how 
to manipulate the boundary region one may reduce drag 
and improve the performance of technological devices. 
The physicist is interested in the small scale structure of 
turbulence away from any boundary, where the action of 
fluid mechanics effectively homogenizes the flow charac- 
teristics and where universal phenomena may be sought. 



In this context "universal" are those phenomena that are 
independent of the nature of the fluid (water, oil, honey 
etc), independent of the mechanism of stirring the flow, 
and independent of the form of the container of the fluid. 
They are inherent features of fluid mechanics as a classi- 
cal field theory. Understanding these universal features 
can gain only marginal technological improvements. But 
this is the theoretical challenge that excites the physicist. 

To many people turbulence research seems orthogo- 
nal to the two main lines of progress in modern physics. 
On the one hand, tremendous effort has been invested 
in understanding the structure of matter, with later de- 
velopments concentrating on ever-diminishing scales of 
constituent particles using the ever-increasing energies of 
particle accelerators. On the other hand astronomy and 
cosmology have exploded with a rich tableau of discov- 
eries at ever-increasing distances from our galaxy. The 
physics of phenomena on the human scale, phenomena 
that are of acute interest to the scientist and layman 
alike, were relegated into a secondary position in the 
course of the development of the first half of 20th cen- 
tury physics. Of course, problems related to the health 
and well-being of humans are deservedly being studied 
in biology and medicine. But physical phenomena that 
can be observed by simply looking out the window are 
considered by many as "non-fundamental" and belong- 
ing to 19th century research. It is the conviction of the 
present writers as well as of a growing number of re- 
searchers that physics on the human scale offers tremen- 
dously rewarding intellectual challenges, some of which 
were at the core of the recent interest in chaotic phenom- 
ena and in the area which is vaguely termed "physics of 
complex systems" . Fluid turbulence, which is the highly 
complex, chaotic and vortical flow that is characteristic 
of all fluids under large stresses, is a paramount example 
of these phenomena that are immensely challenging to 
the physicist and the mathematician alike. The aim of 
this paper is to explain why this problem is exciting, why 
it is difficult, and what are the possible routes that one 
can traverse in finding the solution. The point of view 
described here is that of the physicist whose interest is 
biased in favour of universal phenomena. 



II. SOME HISTORY 

The mathematical history of fluid mechanics begins 
with Leonhard Euler who was invited by Frederick the 
Great to Potsdam in 1741. According to a popular story 
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(which we could not corroborate) one of his tasks was to 
engineer a water fountain. As a true theorist, he began 
by trying to understand the laws of motion of fluids. In 
1755 he wrote Newton's laws for a fluid which in modern 
notation density reads (for the case of constant density) 



du(r, t) 
dt 



u(r, t) ■ Vu{r, t) = -Vp(r, t) 



(1) 



Here u(r, t) andp(r, t) are the fluid velocity and pressure 
at the spatial point r at time t. The LHS of this "Euler 
equation" for u(r, t) is just the material time derivative 
of the momentum, and the RHS is the force, which is 
represented as the gradient of the pressure imposed on 
the fluid. In fact, trying to build a fountain on the basis 
of this equation was bound to fail. This equation pre- 
dicts, for a given gradient of pressure, velocities that are 
much higher than anything observed. One missing idea 
was that of the viscous dissipation that is due to the fric- 
tion of one parcel of fluid against neighboring ones. The 
appropriate term was added to (|l|) by Navier in 1827 by 
Stokes in 1845 |0j- The result is known as the "Navier- 
Stokes equations" : 



du(r, t) 



u(r, t) ■ Vu(r, t) = - Vp(r, t) + vV 2 u(r, t) . 



(2) 



Here v is the kinematic viscosity, which is about 1CP 2 
and 0.15 cm 2 /sec for water and air at room temperature 
respectively. Without the term v\7 2 u(r, t) the kinetic 
energy u 2 /2 is conserved; with this term kinetic energy 
is dissipated and turned into heat. The effect of this term 
is to stabilize and control the nonlinear energy conserving 
Euler equation (|l|). 

Straightforward attempts to assess the solutions of this 
equation may still be very non-realistic. For example, we 
could estimate the velocity of water flow in any one of 
the mighty rivers like the Nile or the Volga which drop 
hundreds of meters in a course of about a thousand kilo- 
meters. The typical angle of inclination a is about 1CP 4 
radians, and the typical river depth L is about 10 meters. 
Equating the gravity force ag (g ~ 10 3 cm/sec 2 ) and the 
viscous drag vd 2 ujdz 2 ~ vu/L 2 we find u to be of the or- 
der of 10 7 cm/sec instead of the observed value of about 
10 2 cm/sec. This is of course absurd, perhaps to the re- 
gret of the white water rafting industry. The resolution 
of this discrepancy was suggested by Reynolds || who 
stressed the importance of a dimensionless ratio of the 
nonlinear term to the viscous term in (|2|) . With a veloc- 
ity drop of the order of U on a scale L the nonlinear term 
is estimated as U 2 /L. The viscous term is about vU/L 2 . 
The ratio of the two, known as the Reynolds number Re, 
is UL/v. The magnitude of Re measures how large is the 
nonlinearity compared to the effect of the viscous dissipa- 
tion in a particular fluid flow. For Re<C 1 one can neglect 



the nonlinearity and the solutions of the Navier-Stokes 
equations can be found in close-form in many instances 
0. In natural circumstances Re is very large. For ex- 
ample in the rivers discussed above Re~ 10 7 . Reynolds 
understood that for Re 3> 1 there is no stable stationary 
solution for the equations of motion. The solutions are 
strongly affected by the nonlinearity, and the actual flow 
pattern is complicated, convoluted and vortical. Such 
flows are called turbulent. 

Modern concepts about high Re number turbulence 
started to evolve with Richardson's insightful contribu- 
tions 0] which contained the famous "poem" that para- 
phrased J. Swift: "Big whirls have little whirls that feed 
on their velocity, and little whirls have lesser whirls and 
so on to viscosity -in the molecular sense" . In this way 
Richardson conveyed an image of the creation of turbu- 
lence by large scale forcing, setting up a cascade of energy 
transfers to smaller and smaller scales by the nonlincar- 
ities of fluid motion, until the energy dissipates at small 
scales by viscosity, turning into heat. This picture led 
in time to innumerable "cascade models" that tried to 
capture the statistical physics of turbulence by assuming 
some thing or other about the cascade process. Indeed, 
no one in their right mind is interested in the full solution 
of the turbulent velocity field at all points in space-time. 
The interest is in the statistical properties of the tur- 
bulent flow. Moreover the statistics of the velocity field 
itself is too heavily dependent on the particular boundary 
conditions of the flow. Richardson understood that uni- 
versal properties may be found in the statistics of velocity 
differences Su(ri, r%) = u(r2)—u{r{) across a separation 
R = r"2 — r±. In taking such a difference we subtract the 
non-universal large scale motions (known as the "wind" 
in atmospheric flows). In experiments (see for example 
) it is common to consider one dimensional cuts of 
the velocity field, Sue(R) = Su(r±,r2)-R/R- The interest 
is in the probability distribution function of Sui(R) and 
its moments. These moments are known as the "struc- 
ture functions" S n (R) = (5ui(R) n ) where (. . .) stands for 
a suitably defined ensemble average. For Gaussian statis- 
tics the whole distribution function is governed by the 
second moment S^i?), and there is no information to be 
gained from higher order moments. In contrast, hydro- 
dynamic experiments indicate that turbulent statistics is 
extremely non-Gaussian, and the higher order moments 
contain important new information about the distribu- 
tion functions. 

Possibly the most ingenious attempt to understand the 
statistics of turbulence is due to Kolmogorov who in 1941 
p] proposed the idea of universality (turning the study 
of small scale turbulence from mechanics to fundamen- 
tal physics) based on the notion of the "incrtial range" . 
The idea is that for very large values of Re there is a wide 
separation between the "scale of energy input" L and the 
typical "viscous dissipation scale" n at which viscous fric- 
tion become important and dumps the energy into heat. 
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In the stationary situation, when the statistical charac- 
teristics of the turbulent flow are time independent, the 
rate of energy input at large scales (L) is balanced by the 
rate of energy dissipation at small scales (77), and must be 
also the same as the flux of energy from larger to smaller 
scales (denoted e) as it is measured at any scale R in the 
so-called "inertial" interval rj -C R <C L. Kolmogorov 
proposed that the only relevant parameter in the inertial 
interval is e, and that L and rj are irrelevant for the sta- 
tistical characteristics of motions on the scale of R. This 
assumption means that R is the only available length for 
the development of dimensional analysis. In addition we 
have the dimensional parameters e and the mass density 
of the fluid p. From these three parameters we can form 
combinations p x e y R z such that with a proper choice of 
the exponents x, y, z we form any dimensionality that we 
want. This leads to detailed predictions about the sta- 
tistical physics of turbulence. For example, to predict 
S n (R) we note that the only combination of e and R 
that gives the right dimension for S n is (ei?)™^ 3 . In par- 
ticular for n — 2 this is the famous Kolmogorov "2/3" 
law which in Fourier representation is also known as the 
"-5 /3" law. The idea that one extracts universal proper- 
ties by focusing on statistical quantities can be applied 
also to the correlations of gradients of the velocity field. 
An important example is the rate e(r, t) at which energy 
is dissipated into heat due to viscous damping. This rate 
is roughly i/|Vu(r, t)\ 2 . One is interested in the fluctua- 
tions of the energy dissipation e(r, t) about their mean e, 
e(r, t) = e(r, t) — e, and how these fluctuations are corre- 
lated in space. The answer is given by the often-studied 
correlation function K ee (R) = (e(r + R, i)e(r, t)). If the 
fluctuations at different points were uncorrelated, this 
function would vanish for all R ^ 0. Within the Kol- 
mogorov theory one estimates K €<i (R) ~ v 2 e 4 ' 3 R~ 8 ' 3 , 
which means that the correlation decays as a power, like 
1/i? 8 / 3 . 

Experimental measurements show that Kolmogorov 
was remarkably close to the truth. The major aspect of 
his predictions, i.e. that the statistical quantities depend 
on the length scale R as power laws is corroborated by ex- 
periments. On the other hand, the predicted exponents 
seem not to be exactly realized. For example, the exper- 
imental correlation K ee (R) decays according to a power 
law, K €£ (R) ~ R~^ for 77 <C R <C L, with p having a nu- 
merical value of 0.2— 0.3 instead of 8/3 ||. The structure 
functions also behave as power laws, S n (R) ~ R (n , but 
the numerical values of ( n deviate progressively from n/3 
when n increases Jt],^| . Something fundamental seems to 
be missing. The uninitiated reader might think that the 
numerical value of this exponent or another is not a fun- 
damental issue. However one needs to understand that 
the Kolmogorov theory exhausts the dimensions of the 
statistical quantities under the assumption that e is the 
only relevant parameter. Therefore a deviation in the 



numerical value of an exponent from the prediction of 
dimensional analysis requires the appearance of another 
dimensional parameter. Of course there exists two di- 
mensional parameters, i.e. L and rj which may turn out 
to be relevant. Indeed, experiments indicated that for the 
statistical quantities mentioned above the energy-input 
scale L is indeed relevant and it appears as a normaliza- 
tion scale for the deviations from Kolmogorov's predic- 
tions: S n (R) ~ {eR) n/3 {L/R) Sn where Q n = n/3 - 5 n . 
Such form of scaling which deviates from the predictions 
of dimensional analysis is referred to as "anomalous scal- 
ing" . The realization that the experimental results for 
the structure functions were consistent with L rather 
than rj as the normalization scale developed over a long 
time and involved a large number of experiments; re- 
cently the accuracy of determination of the exponents 
increased appreciably as a result of a clever method of 
data analysis by Benzi, Ciliberto and coworkers Sim- 
ilarly a careful demonstration of the appearance of L in 
the dissipation correlation was achieved by Sreenivasan 
and coworkers || . A direct analysis of scaling exponents 
Q n and /x in a high Reynolds number flow was presented 
by Praskovskii and Oncley, leading to the same conclu- 
sions |ic|| . 

III. TURBULENCE AS A FIELD THEORY 

Theoretical studies of the universal small scale struc- 
ture of turbulence can be classified broadly into two main 
classes. Firstly there is a large body of phenomenologi- 
cal models that by attempting to achieve agreement with 
experiments reached important insights on the nature of 
the cascade or the statistics of the turbulent fields [ pd| . 
In particular there appeared influential ideas, following 
Mandelbrot |l2] ], about the fractal geometry of highly 
turbulent fields which allow scaling properties that are 
sufficiently complicated to include also non-Kolmogorov 
scaling. Parisi and Frisch showed that by introducing 
multifractals one can accomodate the nonlinear depen- 
dence of C« and n |l3| . However these models are not de- 
rived on the basis of the equations of fluid mechanics; one 
is always left with uncertainties about the validity or rel- 
evance of these models. The second class of approaches is 
based on the equations of fluid mechanics. Typically one 
acknowledges the fact that fluid mechanics is a (classical) 
field theory and resorts to field theoretic methods in order 
to compute statistical quantities. Even though there had 
been a continuous effort during almost 50 years in this 
direction, the analytic derivation of the scaling laws for 
K ee (R) and S n (R) from the Navier-Stokes equations (||) 
and the calculation of the numerical value of the scaling 
exponents p, and £ n have been among the most elusive 
goals of theoretical research. Why did it turn out to be 
so difficult? 
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To understand the difficulties, we need to elaborate a 
little on the nature of the field theoretic approach. Sup- 
pose that we want to calculate the average response of a 
turbulent fluid at some point ro to forcing at point ri. 
The field theoretic approach allows us to consider this 
response as an infinite sum of all the following processes: 
firstly there is the direct response at point ro due to the 
forcing at ri. This response is instantaneous if we as- 
sume that the fluid is incompressible (and therefore the 
speed of sound is infinite). Then there is the process of 
forcing at ri with a response at an intermediate point r 2 , 
which then acts as a forcing for the response at ro. This 
intermediate process can take time, and we need to in- 
tegrate over all the possible positions of point r 2 and all 
times. This is the second-order term in perturbation the- 
ory. Then we can force at ri , the response at r 2 acting as 
a forcing for r 3 and the response at r 3 forces a response 
at ro . We need to integrate over all possible intermediate 
positions r 2 and r 3 and all the intermediate times. This 
is the third-order term in perturbation theory. And so 
on. The actual response is the infinite sum of all these 
contributions. In applying this field theoretical method 
one encounters three main difficulties: 

(A) The theory has no small parameter. The usual 
procedure is to develop the theory perturbatively around 
the linear part of the equation of motion. In other words, 
the zeroth order solution of Eq. (^) is obtained by discard- 
ing the terms which are quadratic in the velocity field. 
The expansion parameter is then obtained from the ratio 
of the quadratic to the linear terms; this ratio is of the- 
order of Reynolds number Re which was defined above. 
Since we are interested in Re3> 1, naive perturbation ex- 
pansions are badly divergent. In other words the contri- 
bution of the various processes described above increases 
as (Re)™ with the number n of intermediate points in 
space-time. 

(B) The theory exhibits two types of nonlinear inter- 
actions. Both are hidden in the nonlinear term u ■ Vm in 
Eq. (En. The larger of the two is known to any person 
who watched how a small floating object is entrained in 
the eddies of a river and swept along a complicated path 
with the turbulent flow. In a similar way any fluctua- 
tion of small scale is swept along by all the larger eddies. 
Physically this sweeping couples any given scale of mo- 
tion to all the larger scales. Unfortunately the largest 
scales contain most of the energy of the flow; these large 
scale motions are what is experienced as gusts of wind 
in the atmosphere or the swell in the ocean. In the per- 
turbation theory for S n (R) one has the consequences of 
the sweeping effect from all the scale larger than R, with 
the main contribution coming from the largest, most in- 
tensive gusts on the scale of L. As a result these con- 
tributions diverge when L — > 00. In the theoretical jar- 
gon this is known as "infrared divergences" . Such diver- 
gences are common in other field theories, with the best 
known example being quantum electrodynamics. In that 



theory the divergences are of similar strength in higher 
order terms in the series, and they can be removed by 
introducing finite constants to the theory, like the charge 
and the mass of the electron. In the hydrodynamic the- 
ory the divergences become stronger with the order of 
the contribution, and to eliminate them in this manner 
one needs an infinite number of constants. In the jargon 
such a theory is called "not renormalizable" . However, 
sweeping is just a kinematic effect that does not lead to 
energy redistribution between scales, and one may hope 
that if the effect of sweeping is taken care of in a consis- 
tent fashion a renormalizable theory might emerge. This 
redistribution of energy results from the second type of 
interaction, that stems from the shear and torsion effects 
that are sizable only if they couple fluid motions of com- 
parable scales. The second type of nonlinearity is smaller 
in size but crucial in consequence, and it may certainly 
lead to a scale-invariant theory. 

(C) Nonlocality of interaction in r space. One recog- 
nizes that the gradient of the pressure is dimensionally 
the same as (u ■ V)ii, and the fluctuations in the pres- 
sure are quadratic in the fluctuations of the velocity. This 
means that the pressure term is also nonlinear in the ve- 
locity. However, the pressure at any given point is de- 
termined by the velocity field everywhere. Theoretically 
one sees this effect by taking the divergence of Eq.(||). 
This leads to the equation V 2 p = V • [(u ■ V)ti]. The 
inversion of the Laplacian operator involves an integral 
over all space. Physically this stems from the fact that in 
the incompressible limit of the Navier-Stokes equations 
sound speed is infinite and velocity fluctuations in all 
distant points are instantaneously coupled. 

Indeed, these difficulties seemed to complicate the ap- 
plication of field theoretic methods to such a degree that 
a wide-spread feeling appeared to the effect that it is im- 
possible to gain valuable insight into the universal prop- 
erties of turbulence along these lines, even though they 
proved so fruitful in other field theories. The present au- 
thors (as well as other researchers starting with Kraich- 
nan [] and recently Migdal JlJ] , Polyakov jl5| , Eyink 
etc.) think differently, and in the rest of this paper we 
will explain why. 

The first task of a successful theory of turbulence is 
to overcome the existence of the interwoven nonlinear ef- 
fects that were explained in difficulty (B). This is not 
achieved by directly applying a formal field-theoretical 
tool to the Navier-Stokes equations. It does not matter 
whether one uses standard field theoretic perturbation 
theory [[17| , path integral formulation, rcnormalization 
group |]1 8| e-expansion, large A^-limit or one's formal 
method of choice. One needs to take care of the partic- 
ular nature of hydrodynamic turbulence as embodied in 
difficulty (B) first, and then proceed using formal tools. 

The removal of the effects of sweeping is based on 
Richardson's remark that universality in turbulence is 
expected for the statistics of velocity differences across a 



4 



length scale R rather than for the statistics of the veloc- 
ity field itself. The velocity fields are dominated by the 
large scale motions that are not universal since they are 
produced directly by the agent that forces the flow. This 
forcing agent differs in different flow realizations (atmo- 
sphere, wind tunnels, channel flow etc.). Richardson's in- 
sight was developed by Kraichnan who attempted to cast 
the field theoretic approach in terms of Lagrangian paths, 
meaning a description of the fluid flow which follows the 
paths of every individual fluid particle. Such a descrip- 
tion automatically removes the large scale contributions 
p0| . Kraichnan's approach was fundamentally correct, 
and gave rise to important and influential insights in the 
description of turbulence, but did not provide a conve- 
nient technical way to consider all the orders of pertur- 
bation theory. The theory does not provide transparent 
rules how to consider an arbitrarily high term in the per- 
turbation theory. Only low order truncations were con- 
sidered. 

A way to overcome difficulty (B) was suggested by 
Bclinicher and L'vov J2f| who introduced a novel trans- 
formation that allowed on one hand the elimination of 
the sweeping that leads to infrared divergences, and on 
the other hand allows the development of simple rules 
for writing down any arbitrary order in the perturbation 
theory for the statistical quantities. The essential idea 
in this transformation is the use of a coordinate frame 
in which velocities are measured relative to the veloc- 
ity of one fluid particle. The use of this transformation 
allowed the examination of the structure functions of ve- 
locity differences S n (R) to all orders in perturbation the- 
ory. Of course, difficulty (A) remains; the perturbation 
series still diverges rapidly for large values of Re, but now 
standard field theoretic methods can be used to reformu- 
late the perturbation expansion such that the viscosity is 
changed by an effective "eddie viscosity" . The theoreti- 
cal tool that achieves this exchange is known in quantum 
field theory as the Dyson line resummation fl23| . The 
result of this procedure is that the effective expansion 
parameter is no longer Re but an expansion parameter 
of the order of unity. Of course, such a perturbation se- 
ries may still diverge as a whole. Notwithstanding it is 
crucial to examine first the order-by-order properties of 
series of this type. 

Such an examination leads to a major surprise: ev- 
ery term in this perturbation theory remains finite when 
the energy-input scale L goes to oo and the viscous- 
dissipation scale r] goes to 23 . The meaning of this 



is that the perturbative theory for S n does not indicate 
the existence of any typical length-scale. Such a length 
is needed in order to represent deviations in the scaling 
exponents from the predictions of Kolmogorov's dimen- 
sional analysis in which both scales L and 77 are assumed 
irrelevant. In other areas of theoretical physics in which 
anomalous scaling has been found it is common that al- 
ready the perturbative series indicates this phenomenon. 



In many cases this is seen in the appearance of logarith- 
mic divergences that must be tamed by truncating the 
integrals at some renormalization length. Hydrodynamic 
turbulence seems at this point different. The nonlinear 
Bclinicher-L'vov transformation changes the underlying 
linear theory such that the resulting perturbative scheme 
for the structure functions is finite order by order |^l|,^3| . 
The physical meaning of this result is that as much as can 
be seen from this perturbative series the main effects on 
the statistical quantities for velocity differences across a 
scale R come from activities of scales comparable to R. 
This is the perturbative justification of the Richardson- 
Kolmogorov cascade picture in which widely separated 
scales do not interact. 

Consequently the main question still remains: how 
does a renormalization scale appear in the statistical the- 
ory of turbulence? 

It turns out that there are two different mechanisms 
that furnish a renormalization scale, and that finally both 
L and rj appear in the theory. The viscous scale r\ ap- 
pears via a rather standard mechanism that can be seen 
in perturbation theory as logarithmic divergences, but 
in order to see it one needs to consider the statistics of 
gradient fields rather than the velocity differences them- 
selves |2^j2^]. For example, considering the perturbative 
series for K ee (R), which is the correlation function of the 
rate of energy dissipation ^|Vm| 2 , leads immediately to 
the discovery of logarithmic ultraviolet divergences in ev- 
ery order of the perturbation theory. These divergences 
are controlled by an ultraviolet cutoff scale which is iden- 
tified as the viscous-dissipation scale rj acting here as the 
renormalization scale. The summation of the infinite se- 
ries results in a factor (R/i]) 2A with some anomalous 
exponent A which is, generally speaking, of the order of 
unity. The appearance of such a factor means that the 
actual correlation of two i?-separated dissipation fields 
is much larger, when R is much larger than 77, than the 
naive prediction of dimensional analysis. The physical 
explanation of this renormalization |^,^8| is the effect of 
the multi-step interaction of two i?-separated small ed- 
dies of scale T] with a large eddy of scale R via an infinite 
set of eddies of intermediate scales. The net result on the 
scaling exponent is that the exponent /1 changes from 8/3 
as expected in the Kolmogorov theory to 8/3 — 2 A. 

At this point it is important to understand what is the 
numerical value of the anomalous exponent A. In |24| l 
there was found an exact sum that forces a relation be- 
tween the numerical value of A and the numerical value 
of the exponent (2 of 5*2 (R) , A = 2 — (2 ■ Such a relation 
between different exponents is known in the jargon as a 
"scaling relation" or a "bridge relation" . Physically this 
relation is a consequence of the existence of a universal 
nonequilibrium stationary state that supports an energy 
flux from large to small scales |24|,|25|]. The scaling rela- 
tion for A has far-reaching implications for the theory of 
the structure functions. It was explained that with this 
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value of A the series for the structure functions S n (R) 
diverge when the energy-input scale L approaches oo as 
powers of L, like (L/R) Sn . The anomalous exponents S n 
are the deviations of the exponents of S n (R) from their 
Kolmogorov value. This is a very delicate and impor- 
tant point, and we therefore expand on it. Think about 
the series representation of S n (R) in terms of lower order 
quantities, and imagine that one succeeded to resum it 
into an operator equation for S n (R). Typically such a 
resummed equation may look like [1 — 0]S n (R)] =RHS, 
where O is some integro-differential operator which is 
not small compared to unity. If we expand this equa- 
tion in powers of O around the RHS we regain the infi- 
nite perturbative series that we started with. However, 
now we realize that the equation possesses also homoge- 
neous solutions, solutions of [1 — 0]S n (R)] = which are 
inherently nonperturbative since they can no longer be 
expanded around a RHS. These homogeneous solutions 
may be much larger than the inhomogeneous perturba- 
tive solutions. Of course, homogeneous solutions must 
be matched with the boundary conditions at R = L, and 
this is the way that the energy input scale L appears in 
the theory. This is particularly important when the ho- 
mogeneous solution diverge in size when L — > oo as is 
indeed the case for the problem at hand. 

The next step in the theoretical development is to un- 
derstand how to compute the anomalous exponents S n . 
The divergence of the perturbation theory for S n (R) with 
L — > oo forces us to seek a nonperturbative handle on the 
theory. One finds this in the idea that there exists always 
a global balance between energy input and dissipation, 
which may be turned into a nonperturbative constraint 
on each n-th order structure function [25| . Using the 
Navier-Stokes equations (|J) one derives the set of equa- 
tions of motion 
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where D n and J n stem from the nonlinear and the vis- 
cous terms in (|^) respectively. To understand the physi- 
cal meaning of this equation note that S*2(i?) is precisely 
the mean kinetic energy of motions of size R. The term 
D2(R) whose meaning is the rate of energy flux through 
the scale R is known exactly: D2(R) — dS^i?) / dR. The 
term vJ2{R) is precisely the rate of energy dissipation 
due to viscous effects. The higher order equation for 
n > 2 are direct generalizations of this to higher order 
moments. In the stationary state the time derivative van- 
ishes and one has the balance equation D n (R) — v,J n (R) . 
For n = 2 it reflects the balance between energy flux and 
energy dissipation. The evaluation of D n (R) for n > 2 
requires dealing with the difficulty (C) of the nonlocality 
of the interaction, but it does not pose conceptual diffi- 
culties. It was shown |^5| that D n (R) is of the order of 
dS n +i/dR. On the other hand, the evaluation of J n (R) 
raises a number of very interesting issues whose resolu- 



tion lies at the heart of the universal scaling properties 
of turbulence. Presently not all of these issues have been 
resolved, and we briefly mention here some ground on 
which progress has been made by the present authors. 

From the derivation of Eq.(||) one finds that J n (R) con- 
sists of a correlation of V 2 it with n — 2 velocity differences 
across a scale R = |ri — r 2 )|: (V 2 u(ri)[6ug(ri, r 2 )]"~ 2 ). 
The question is how to evaluate such a quantity in terms 
of the usual structure functions S n (R). Recall that a gra- 
dient of a field is the difference in the field values at two 
points divided by the separation when the latter goes to 
zero. In going to zero one necessarily crosses the dissipa- 
tive scale. To understand what happens in this process 
one needs first to introduce many-point correlation func- 
tions of a product of n velocity differences: 

F n (ro\ri > ...r n ) = (Sui(r ,ri)...6ue(ro,r n )) . (4) 

Next we need to formulate rules for the evaluation of such 
correlation functions of velocity differences when some of 
the coordinates get very close to each other. For example, 
a gradient d/dr a can be formed from the limit v\ — > r$ 
when we divide by ri, a — ro, a . These rules are known 
in the theoretical jargon as "fusion rules". The fusion 
rules for hydrodynamic turbulence were presented in p9| . 
They show that when p coordinates in F n are separated 
by a small distance r, and the remaining n—p coordinates 
are separated by a large distance R, then the scaling 
dependence on r is like that of S p (r), i.e. . This is true 
until r crosses the dissipative scale. Assuming that below 
the viscous-dissipation scale rj derivatives exist and the 
fields are smooth, one can estimate gradients at the end 
ofthe smooth range by dividing differences across n by n. 
The question is, what is the appropriate cross-over scale 
to smooth behaviour? Is there just one cross-over scale 
77, or is there a multiplicity of such scales, depending on 
the function one is studying? For example, when does 
the above n-point correlator become differentiable as a 
function of r when p of its coordinates approach Tq7 Is 
that typical scale the same as the one exhibited by S p (r) 
itself, or does it depend on p and n and on the remaining 
distances of the remaining n — p coordinates that are still 
far away from ro ? 

The answer is that there is a multiplicity of cross-over 
scales. For the n-point correlator discussed above we de- 
note the dissipative scale as rj{p,n,R), and it depends on 
each of its arguments |29|,^6| . In particular it depends on 
the inertial range variables R and this dependence must 
be known when one attempts to determine the scaling ex- 
ponents ( n ofthe structure functions. In brief, this line of 
thought leads to a set of non-trivial scaling relations. For 
example we confirm the phenomenologically conjectured 
p"l| "bridge relation" fx — 2 — (in close agreement with 
the experimental values) and predict that although the 
C„ are not Kolmogorov, they are nevertheless linear in n 
for large n. 
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IV. SUMMARY 



physics and engineering will be the challenge of the 21st 
century. 



It appears that there are four conceptual steps in 
the construction of a theory of the universal anomalous 
statistics of turbulence on the basis of the Navier-Stokes 
equations. First one needs to take care of the sweeping 
interactions that mask the scale invariant theory |^l|,^3| . 
After doing so the perturbation expansion converges or- 
der by order, and the Kolmogorov scaling of the velocity 
structure functions is found as a perturbative solution. 
Secondly one understands the appearance of the viscous- 
dissipation scale 77 as the natural normalization scale in 
the theory of the correlation functions of the gradient 
fields [ p7|j24| ]. This step is similar to critical phenomena 
and it leads to a similarly rich theory of anomalous be- 
haviour of the gradient fields. Only the tip of the iceberg 
was considered above. In fact when one considers the 
correlations of tensor fields which are constructed from 
du a /drp (rather than the scalar field e) one finds that ev- 
ery field with a different transformation property under 
the rotation of the coordinates has its own independent 
scaling exponent which is the analog of A above [j30[ . 
The third step is the understanding of the divergence 
of the diagrammatic series for the structure functions as 
a whole (^5|. This sheds light on the emergence of the 
energy-input scale L as a normalization length in the 
theory of turbulence. This means that the Kolmogorov 
basic assertion that there is no typical scale in the ex- 
pressions for statistical quantities involving correlations 
across a scale R when r\ <g; R <C L is doubly wrong. In 
general both lengths appear in dimensionless combina- 
tions and change the exponents from the predictions of 
dimensional analysis. Examples of correlation functions 
in which both normalization scales L and 77 appear simul- 
taneously were given explicitly |30|. Last but not least 
is the formulation of the fusion rules and the exposition 
of the multiplicity of the dissipative scales which should 
eventually result in a satisfactory description of all the 
scaling properties p9 |. 

The road ahead is not fully charted, but it seems that 
some of the conceptual difficulties have been surmounted. 
We believe that the crucial building blocks of the the- 
ory are now available, and they begin to delineate the 
structure of the theory. We hope that the remaining 4 
years of this century will suffice to achieve a proper un- 
derstanding of the anomalous scaling exponents in tur- 
bulence. Considerable work, however, is still needed in 
order to fully clarify many aspects of the problem, and 
most of them are as exciting and important as the scaling 
properties. There are universal aspects that go beyond 
exponents, such as distribution functions and the eddy 
viscosity, and there are important non-universal aspects 
like the role of inhomogeneities, the effect of boundaries 
and so on. Progress on these issues will bring the theory 
closer to the concern of the engineers. The marriage of 
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